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Abstract 

We present results of a perturbative QCD analysis of deep inelastic measurements of both the 
deuteron and proton structure functions. We evaluate the theoretical uncertainty associated to 
nuclear effects in the deuteron, and we extract simultaneously the isospin depedendence of: i) 
the higher twists terms; ii) the ratio of the longitudinal to transverse cross sections, R = ai^jaT'^ 
in) the ratio of the neutron to proton structure functions, -Fl^/Fl* . The extraction of the latter, 
in particular, has been at the center of an intense debate. Its accurate determination is crucial 
both theoretically and for the interpretation of the more precise neutrino experiments including 
the newly planned high intensity 50 GeV proton synchrotron. 
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I. INTRODUCTION 



Deep Inelastic Scattering (DIS) experiments provide the most accurate measurements of 
the strong couphng constant, as, at intermediate scales. They are also the main source of 
information on the Parton Distribution Functions (PDFs) in the proton and neutron. The 
precision with which both of these quantities are known reflects directly into the precision of 
calculations of the cross sections for all other hard scattering processes. An accurate deter- 
mination of these, in turn, plays a key role both in the extraction of possible contributions 
of new physics at new collider energies, and in the interpretation of the forthcoming high 
precision experiments using neutrino beams 0. 

Both as and PDFs are not directly observable and they need to be extracted from the 
DIS data according to some procedure. A number of uncertainties affect the analysis, re- 
lated to both the perturbative QCD (pQCD) series - inclusion of higher orders, threshold 
resummation effects - and to corrections that are non-perturbative in nature - Target Mass 
Corrections (TMC), dynamical power corrections, and nuclear effects in the case of the 
neutron structure function. 

It is therefore mandatory to be able to control the size of these uncertainties by intro- 
ducing a systematic, well tested, method of extraction in which possible ambiguities can be 
properly gauged. While analyses along these lines exist for the proton structure function 
([3] in the DIS region and jsj in the "few GeV" or resonance region), an accurate and com- 
plete treatment of the neutron structure function is still lacking. This paper is devoted to 
the application of a newly developed method to determine the isospin dependence of the 
nucleon Structure Functions (SFs). With the analysis presented here we hope to contribute 
to the interpretation of both recent data and new experiments, by providing a quantitative 
measure of the space in which pQCD based corrections and nuclear effects can be wiggled. 

For our analysis we use the extensive DIS measurements that exist for both proton and 
deuteron targets in a wide range of kinematics with the exception for the very large x region 
where fewer data sets exist, mostly at low values of the final state invariant mass, W"^, in 
the region of nucleon resonances (more experiments are however beingplanned that will 
cover the large x DIS region in forthcoming programs at Jefferson Lab and at neutrino 
facilities P). 

In QCD, different contributions to the DIS structure functions can be written using the 
Operator Product Expansion (OPE), by ordering them according to their twist, r (r = 
dimension — spin) j^. The Leading Twist (LT) contribution (with r = 2 in DIS) is di- 
rectly related to the single particle properties of quarks and gluons inside the nucleon - the 
PDFs. The Higher Twist (HT) components (r = 4,6, ... in unpolarized DIS) involve inter- 
actions between quarks and gluons in the nucleon and they are suppressed by terms of order 
l/Q^, 1/Q^..., respectively. In phenomenological studies, the PDFs are extracted from QCD 
global fits. Accurate extractions use data with sufficiently high and invariant mass W'^, 
where both target mass and HT corrections are expected to be very small. QCD fits can 
now be performed to order a| (Next-to-Next-to-Leading Order, NNLO approximation). If 
the data encompass a large range in Q^, higher order corrections as well as HT effects need 
to be taken into account simultaneously. 

In a recent series of papers 0] a proper choice of the statistical estimator allows one to 
propagate all experimental error into the uncertainties in the PDFs. Because of the statis- 
tical efficiency of this new estimator, the overall systematic error on the PDFs is sensibly 
reduced with respect to previous analyses based on simplified estimators. With a better de- 
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termination of systematic experimental errors in hand, one can address in detail the sources 
of theoretical errors. Theoretical uncertainties are in principle an elusive concept as by 
definition they refer to quantities that have yet to be calculated. Uncertainties/ambiguities 
of this type and inherent to the pQCD analyses are due to: i) impact of the higher order 
QCD corrections; ii) the HT terms; in) target mass corrections (TMC); iv) heavy quarks 
masses and threshold values; v) the form of the initial PDFs. These questions were ad- 
dressed in detail in the analysis of 0, where the DIS cross sections were fitted by including 
both the LT terms calculated from the PDFs evolved to the NNLO QCD, and the twist-4 
terms, evaluated separately for the proton and neutron structure functions F2 and Fl- The 
fact that the cross sections data were fitted instead of the data on F2, allowed for a better 
determination of Fl, which in j3] was obtained iteratively. In summary, the analysis of 
shows that theoretical uncertainties from the pQCD series are under control, and that, due 
to the new estimator, all extracted quantities can be determined with smaller errors than in 
previous analyses. 

An additional uncertainty is however present in analyses of the isospin dependence of 
both PDFs and the HT terms, in that the neutron structure functions have to be extracted 
from nuclear data. The main thrust of our analysis has been to make a thorough assessment 
of the impact of nuclear corrections on both the LT and HT terms. Here, in particular, 
we focus on the isospin dependence of the HT terms. Detailed results on both the ratio 
F2/F2 and on R = aL/c^T will be presented in a forthcoming paper In our analysis 
we use the deuteron data where uncertainties are expected to be in better control. We 
address uncertainties arising from: i) Different models of nuclear effects - we highlight in 
particular the differences with using the extrapolation 01 of the nuclear density model of 
the EMC effect jsj; ii) Different deuteron wave functions derived from currently available 
NN potentials, giving rise to different amounts of high momentum components; Hi) The 
interplay between nucleon off-shellness and TMC in nuclei. 

In addition to its practical purpose, a quantitative determination of the isospin depen- 
dence of the HTs contributions is of theoretical interest in understanding the nature of power 
corrections. On one side, infrared renormalons (IRR) have been suggested as a method for 
estimating the contribution of power corrections to the cross sections for a number of hard 
processes (see j9| and references therein). Based on this hypothesis the calculations in 
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have predicted the x dependence of the coefficients of the HT terms for both F^ and the 
valence and singlet contributions to F2. On the other hand, some models exist that predict a 
sizeable isospin dependence of the HT terms. This has been suggested for instance in models 
that interpolate between partonic and non-partonic degrees of freedom at low as in 3 



(a smaller effect seems to occur, however, in the predictions for the HT isospin dependence 
in 01). A large effect of about a factor two for the ratio of the neutron to proton HT terms. 



was also predicted based on quark counting estimates in Ref. [14 



A thorough analysis of the isospin dependence of HTs might therefore help disentangling 
predictions of different models, by investigating, for instance, whether the onset of "flavor- 
blind" power corrections from IRR occurs at a different scale than for other dynamical 
ones. 

In what follows, we present our results for each quantity along with a discussion of the 
nature of the nuclear effects and the extraction method. In Section |H] we present the general 
formalism and deflnitions. In Section IIIII we discuss the contribution of nuclear effects. In 
Section IIVI we outline the extraction method and we present our results. In Section |3 we 
discuss some phenomenological applications of our analysis. Finally, we draw our conclusions 
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in Section IVTl 



II. DIS FORMALISM 

The inclusive DIS section of unpolarized charged leptons off an unpolarized nucleon or 
nuclear target is fully determined by the spin-averaged electromagnetic tensor of the target, 
which can be parameterized in terms of two invariant structure functions Fi and F2 (we 
ignore a small contribution due to neutral current) 

W,Ap,q) = / d'zexp{tq-z){p,s\ [j;'^{z), J^iO)] \p,s) 

s 

= -9,uF, + ^F„ (1) 
p-q 

where J^™ is electromagnetic current, p and q are the target momentum and four-momentum 
transfer, respectively. In order to simplify notations, we denote 

9^lu — g^iu ^ 
~ P-q /ON 

= P/. - ^g,.. (3) 

The normalization of states adopted here is {p\p') = 2j)o(2vr)^5(p — p') for both the bosons 
and fermions. With this normalization the structure functions F12 are dimensionless. They 
depend on two invariant variables, namely the Bjorken variable x = Q'^/2p ■ q and the 
four- momentum transfer squared = —q^. 

The differential cross section in terms of the structure functions and standard variables 
X and y = p ■ q/p ■ k, where k is the incoming lepton four-momentum, reads 

d^a ina' 



dxdy Q'^xy 



{Mxyf\ ^ 1 2 / 



(4) 



where a is the electromagnetic coupling constant, M is the nucleon mass and = 2xyp- k , 
and Ft = 2x^1. 0. 

In QCD the contributions from different quark-gluon operators to the electromagnetic 
tensor and to DIS structure functions are ordered according to their twist, leading to the 
expansion in inverse powers of Q^: 

FtAx, Q') = Ft^J (x, Q') + + O (1/g^) . (5) 

This expansion applies to both the proton and neutron structure functions, we have sup- 
pressed the indices for simplicity. The first term is the LT contribution and H12 are the HT 
- twist-4 - contributions. Furthermore, if a finite mass for the nucleon target is considered, 
new terms arise in Eq.© that mix operators of different spin, leading to additional power 
terms of kinematical origin - the so-called target mass corrections. In the approximation 
that x^M^/Q^ is small, the TMC series can be absorbed in the leading twist term [l5l |. 

The separation of LT, TMC and dynamical HT from the data, is not straightforward as 
witnessed by the number of studies dedicated to it since the initial formulation of the problem 
m the 1970's HI, Eel. In what follows we define the LT, TMC, and HT, contributions to the 



structure functions Ft and F2. The formalism for nuclear DIS is discussed in Section ITTTl 



4 



A. Leading Twist 



The LT part of the structure functions is related to the PDFs, pi{x,Q) - the index i 
refers to the different types of quarks and antiquarks, and to the gluon distribution - via a 
convolution with perturbatively calculable coefficient functions C^2- 

/■I 

F^l{x,Q') = J2 / -C'T,2{z,c.sm)M^/^.Q')- (6) 

The dependence of the PDFs is predicted by the well-known DGLAP evolution equa- 
tions Q: 

E f^P.,{z,as{t))p,{x/z,t), (7) 

j=q,q,g ^ 

where t = Q^, and as is the strong coupling constant, and Pij are the splitting functions. 

The coefficient functions have been calculated to NNLO The splitting functions 

are known to NLO, and only a limited set of Mellin moments are evaluated to NNLO j2ol |. 
Although, estimates of the full x-dependence of the splitting functions in NNLO approxi- 
mation are available in our analysis we use the MS NLO QCD approximation with the 
renormalization/factorization scales chosen equal to Q. The NNLO variant of our fit is used 
to estimate the uncertainty due to higher orders. Large x resummation effects, arising from 
terms of the type [asiQ"^) ln(l — z)]^'^ at in the coefficient functions are present in 

principle. They have been shown to be comparable in size to NNLO corrections in Ref. 
and to generate a further negative correction to the HT coefficient of the proton SF. 
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B. Target Mass Corrections 

We follow the method of Ref. where it was shown that the TMC series can be summed 
up, leading to the modification of the LT term. Therefore, Eq.(0) remains valid with the LT 
terms replaced by 

F-™%, Q^) = Q') + 2^^ Q% (8a) 

F^™^ix, = Q') + Q% (8b) 

where 7 = (1 + ix'^M'^ /Q'^Y^^ and C, = 2x/ (1 + 7) is the Nachtmann variable 0. 

It must be noted however, that the derivation of 0] was given in the zeroth order in as, 
assuming that the target quarks are on-shell. Both, higher order as corrections and quark 
off-shell effects modify Eqs.®- It was argued in [23[ that off-shell effects lead to M^/Q^ 
terms which are not incorporated in Eqs.(jHl). addition, target mass corrections should be 
applied also to the HT terms in the higher order terms in the expansion (jSJ. For this reason 
we do not consider l/Q"^ terms in the TMC formula, which are small for the kinematical 
range considered. Finally, TMC corrections for an off-shell target, i.e. when 7^ M^, 
should be treated as part of the nuclear effects and will be discussed therefore in Section III. 
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C. Higher Twists 



The extraction of higher twist terms from the data is a longstanding problem, as recog- 
nized from the very first developments of a pQCD phenomenology (0, 2J|)- HTs have been 
hard to pin down for a number of reasons. First of all a connection with partonic interpreta- 
tions cannot be established on a one by one basis, differently from the LT components that 
are directly related to the PDFs. In fact the HT terms are formally written within OPE as 
the product of coefficient functions and hadronic matrix elements of composite local opera- 
tors. Not all of the matrix elements are independent, but a minimal basis can be selected 
after relating them through the equations of motion j2^. Nevertheless, even in a minimal 
basis, the number of independent reduced matrix elements is much larger than the number 
of observables - e.g. the moments of the structure functions. In unpolarized scattering 
one can single out formally the four quarks and two quarks-two gluons types of operators, 
corresponding in a partonic language to quark-quark and quark-gluon correlations, respec- 
tively. The determination of the relative scales of these contributions is at present model 
dependent. However, it was shown in Ref. [2^ that a simultaneous analysis of both F2 and 
Fl better constrained the evaluation of the quark-gluon term, that was required to be large 
and negative in order to fit both sets of data. 

A practical difficulty is also in that theoretical estimates - from the simple estimate of 
the increase of the number of operators with respect to twist-2 ji^l, to more sophisticated 
IRR calculations [13, 11| - predict that HTs are most important at low W"^ ~ — x). 



In this limit, with > 4 GeV^, thus avoiding the resonance region, it can be easily shown 
that the logarithmic dependencies characterizing pQCD evolution to a given order mimic 
the dependence of the twist-4 terms. A correlation between the pQCD parameters and 
the HT coefficients arises that has lead the authors of jl^l to conclude that for the structure 
function F3, the NNLO term and the HT corrections are, within the precision of current 
data, undistinguishable. A similar analysis was subsequently applied in Ref.01 to F2. It 
lead the authors to the conclusion that HTs are highly reduced with respect to previous 
determinations, even inside the resonance region. 

To summarize, the results of a number of analyses of the HT contributions in DIS are 
still not conclusive. A joint search using a combination of predictions from hadronic models 
and an accurate experimental extraction seems to be the most promising avenue. It is in 
this spirit that in this paper we address yet another aspect of the phenomenology of HTs, 
namely their isospin dependence. An isospin dependence of the HT terms resides entirely 
in their non-perturbative structure. In IRR models in fact the coefficients of the HT terms 
are predicted to have no target dependence, provided the radiative corrections to these 
terms are factored out by assuming for instance, that they behave similarly to their LT 
correspondent. An assessment of the magnitude of the isospin dependence of the HT terms 
provides therefore a handle on understanding their nature and more precisely the extent to 
which they can be described by models. In particular, possible scenarios about the large x 
structure of the proton can be investigated, envisaging coherent scattering from multiquark 
composites carrying increasing momentum fraction at x — > 1 3 • 



III. NUCLEAR EFFECTS 

Since nuclear data are used in the QCD analysis, nuclear effects must be taken into 
account. In our analysis we use deuterium data in order to study the isospin effect in PDFs 
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and HT terms. In this section we briefly discuss our method to correct for nuclear effects. 
More details can be found elsewhere lal. 



A. Fermi motion and binding effects 

For large x > 0.1, away from nuclear shadowing region, nuclear DIS of leptons off nuclear 
targets can be viewed as incoherent scattering off bound nucleons HHHIslllllal. The 



DIS cross section is given by the imaginary part of the virtual photon Compton amplitude 
in the forward direction. In incoherent scattering approximation, the nuclear Compton am- 
plitude is taken in impulse approximation by disregarding both initial state interactions and 
interactions between the struck quark and the nuclear debris. The corresponding Feynman 
diagram is shown in Fig^ Calculating this diagram and projecting the structure functions 
from the imaginary part of the Compton amplitude we derive the relations between the 
nuclear and the nucleon structure functions. For the deuteron we have 

-FFM') = + f ) (9a) 

I'F.^M') = [ ^,\^D{p)fU + ^-^) (9b) 



(27r)3 ' V M 

'^^^Vi^(x',Q^;/), 



where F^^^ = + F^j.)/2. The deuteron wave function \1/d(p) squared describes the 

probability to find the bound proton (or neutron) with momentum p, x' = Q'^/2p-q is the 
Bjorken variable of the bound nucleon with the four-momentum p, which is given by the 
difference of the target four-momentum and the four-momentum of the spectator nucleon. 
Eqs.® are written for the target rest frame and the z axis is chosen such that the momentum 
transfer, q = (go, 0_l, — |q|). In this reference frame p = {Mo — \/p'^ + M'^,p±,Pz), where 
Mo and M are the deuteron and the nucleon mass respectively. The kinematical factors 
in Eqs.(jn)) result from the projection of the structure functions from the hadronic tensor. 
The transverse motion of the bound nucleon in the deuteron rest frame is the reason for the 
appearance of additional terms in the transverse and the longitudinal cross sections, i.e. the 
terms proportional to x''^p\F2 in Eqs.dH)). 



B. OfT-shell effects 

The bound proton and neutron are off-mass-shell and their structure functions differ 
from those of the free proton and neutron. The off-shell nucleon structure functions depend 
on the nucleon virtuality p^ as an additional variable. Therefore, off-shell effects in the 
structure functions are closely related to the target mass corrections. Target mass effects in 
the off-shell nucleon can be of two different kinds. First, similarly to the on-shell nucleon, 
we have to take into account the kinematical target mass dependence due to the finite p^ jQ^ 
ratio. We assume that this effect is described by Eqs.ljHJ, where the nucleon mass squared is 
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replaced by (this leads in turn to the modification of the parameter 7 and the variable ^ 
in the off-shell region). Furthermore, the dependence on appears already at leading twist 
js^, 113, • In order to estimate the off-shell effect in the LT structure functions, we 
consider the quark distribution in terms of the spectral representation S^, 3^, S^l 



g(x,p2) = y ds y dA:V(s,A;2,a;;/), (10) 
kl^Ax.s) = xlp"--^. (11) 



X 



The integration in Eq.()10|) is taken over the mass spectrum of spectator states s and quark 
virtuality fc^, k'^^^ is the kinematical maximum of k"^ for the given s and p^. The invariant 
spectral density p measures the probability to find in a nucleon with momentum a quark 
with light-cone momentum x and virtuality k"^ and the remnant system in a state with 
invariant mass s. 

We observe from Eq. ljlUj) that the dependence of quark distributions has two primary 
sources: the one in the upper limit of k"^ integration (kinematical off-shell dependence), and 
an explicit p^ dependence of the quark spectral function p (dynamical off-shell dependence). 

The kinematical off-shell effect causes a negative correction to the bound nucleon structure 
functions and produces an enhanced EMC effect, as first noticed in [s^,!!^. However, if only 
the kinematical off-shell effects are taken into account the number of valence quarks in the 
nucleon would change with p^. It can be seen directly from Eq. ljlUj) that the normalization 
of the quark distribution decreases as p^ decreases, provided that the spectral density is 
positively defined. This effect leads to an overall 1 — 2% depletion of valence quarks in 
the deuteron. Furthermore, the magnitude of this effect increases in heavy nuclei, since the 
average shift from the mass shell of the bound nucleon increases. This observation indicates 
that an off-shell effect of dynamical origin must also be present. 

Dynamical off-shell effects can be viewed as a measure of the nucleon's deformation inside 
the nuclear medium. One possible way to evaluate dynamical off-shell effects is to require 
the conservation of the valence quark number in the nucleon also in the off-shell region 

d 

dxgvai(a;;p2) = 0, (12) 



dp^ 

This equation makes it possible to estimate the off-shell effect minimizing the model depen- 
dence j33| . It was shown that Eg. (1121) results in a partial cancellation between kinematical 
and dynamical off-shell effects j33|,|3J]. However, the off-shell effect in the structure functions 
remains an important correction. 



C. Results and comparison with other approaches 

The effect of nuclear corrections is illustrated in Fig|21 The ratio i?^(a;, Q^) = 
F^{x, Q"^) / F2 {x, Q"^) was calculated by Eqs.dHI) using different approximations. The dotted 
curve corresponds to the standard assumption that the bound proton and neutron structure 
functions in Eqs.ijH)) are equal to those of the free nucleon ones. The competition between 
nuclear binding and Fermi motion determines the shape of R2- In particular, the depletion 
of nuclear structure functions at a; < 0.7 is due to the effect of nuclear binding, while the 
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rise of the R2 ratio at large x is due to the effect of the nucleon momentum distribution 
(Fermi motion) [2^. 

The dashed curve corresponds to the results with the target mass corrections. We notice 
that TMC is an important correction at large be directly seen from Eqs.(jH|). This 

correction modifies the shape and the magnitude of the LT structure function at large x. 
This in turn leads to the softening of the ratio i?^ at large x. 

The solid curve stands for the full calculation with TMC and off-shell effect taken into 
account. We observe, that the off-shell effect is most important in the binding region [x 
between 0.3 and 0.7), where it causes a negative correction to the bound nucleon structure 
functions. 

The region of large x corresponds to small masses of produced hadronic states W"^ = 
-I- Q'^{l/x — 1). For instance, the events with Q'^ = 10 GeV^ and x > 0.75 fall into the 
resonance region. For this reason the DIS parametrization of structure functions, which are 
used in computing the nuclear effects with Eq.Q, are questionable at large x. In order to 
avoid the resonance region and elucidate nuclear effects in the DIS regime, we apply a cut 
at W = 1.8 GeV. The ratio i?^ as a function of W was calculated for a few different values 
of Q^. The results are shown in the lower panel of FigI21 

The off-shell effect is much less important in the deuteron than in heavy nuclei. The 
strength of this effect is governed by the average off-shellness of the bound nucleon A = 
(p2 _ j^2^ i]^2 ^ 2{e - T) /M, where T = (p^) /2M and e = (po) - M are the average kinetic 
and separation energy. In order to illustrate the strength of nuclear binding and off-shell 
effects, we have calculated T, e, and A averaged over the nuclear spectral function for a 
number of nuclei. In FigEl these parameters are plotted as a function of the nuclear masss 
number A. It can be seen that A increases by a factor of five when going from the deuteron 
to heavy nuclei. 

A phenomenological model of the EMC effect in the deuteron was given in Ref . js^ . The 
model is based on an extrapolation of SLAC data on the EMC effect for a number of nuclei 
from ^He to ^^^Au, where the key assumption was made that the quantity {F.^/F.^ ~ 1) ^ 
being the nuclear structure function per nucleon - scales with the nuclear density: 

Fi'/F{--1_,, (13) 



Ff/F« - 1 



where pd and pa are the number densities in the deuteron and in a heavy nucleus (Nuclear 
Density Model (NDM) j^). It was also assumed that this ratio is independent of Q^. Using 
Eq.dinD, the ratio F^/F^^ was obtained in Ref. js^ in terms of the experimentally measured 
ratio F2 / F2 by using: 



F^ ^ PA-Pd \Fi 



(14) 



and by averaging the quantity appearing on the r.h.s. over the SLAC nuclear data. The 
values of F^/F,^ were given for x corresponding to the data bins. These results are shown 
in Fig|21by a shaded region. 

We observe that the NDM attempts to extrapolate "density scaling" to the region of light 
nuclei, where the notion of density is ill defined [s^]. The value pd = 0.024 fm~'^ was used in 
Ref. [36|, which was derived using the r.m.s. radius of the deuteron. However, it is not clear 
what volume is occupied by the deuteron and for this reason pd has a large "theoretical" 
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uncertainty. Since the quantity / — 1 is proportional to pd, this theoretical uncertainty 
will directly translate into an uncertainty for the extrapolated ratio F^/F^ . This was not 
given in Ref.js^ and it is likely to be larger than the errors shown in FigO 



IV. QCD ANALYSIS AND FIT 
A. Fitting Procedure 

In our analysis the data on charged lepton DIS off proton and deuterium targets from 
BCDMS, NMC, HI, ZEUS, and SLAG experiments were used Q. The cut > 2.5 GeV^ 
was applied in order to avoid the region where as is large and the higher-order QCD radiative 
corrections can be out of control. The HERA data with > 250 GeV^ were not used in the 
analysis since the impact of those data on the fit is marginal because of large experimental 
errors. The maximum x in the data set is about 0.9 and minimum center-of-mass energy W 
is about 1.8 GeV (see FigE)). The total number of data points (NDP) is 1381 for the proton 
and 998 for the deuteron; /NDP =1.1 for the best fit. 

In our fit we make use of Eq. (@)) with the LT structure functions given in Eq. ^ corrected 
for the target mass effect by Eq.(jH)). The higher twist terms i?2.T of Eq.(0), the parton distri- 
butions, and the value of as were simultaneously fitted to the data. The parton distributions 
were parametrized in the form used in the earlier analyses of Ref . with initial scale of the 
QCD evolution Qq = 9 GeV^. The evolution equations were solved numerically by direct 
integration in x-space. Our procedure is in agreement with the benchmarks introduced in 
Ref. [3^ that require the precision of the solution to be much better than the accuracy of 
the data used in the analysis. 

The implementation of target mass and nuclear corrections in our fit is as follows. The 
proton and neutron structure functions were calculated as 

Fi^- = T,,T {f^'t"^'''^ } + m:T/Q', (15) 

where we have expressed the target mass correction formula given in Eq.® in terms of the 
functional: 72,t- 

The deuteron structure functions were calculated as 

F^j. = ^2,T {-^:^T + -^2"t} 5 (16) 

where J-'2,t are the nuclear smearing functionals corresponding to Eq.(jni), and Fl"'-^ are the 
proton and neutron structure functions of Eq. lfT^ . However, the implementation of this 
approach slows down the numerical calculation because of 4-dimensional integrations in the 
deuteron cross sections in terms of the QCD-evolved PDFs. We are able to reduce the 
calculation time by using approximate expressions for the deuteron structure functions: 



2,T 



2.T 



^piQPM) 
2.T 



j,n{QPM) I 



2.T 



ppiQPM) pn(QPM) 



X 



2,T 



,p{QPM) _^ pniQPM) \ 



2,T 



2,T 



^p{QPM) ^ pn{QPM) 



2,T 



2,T 



2,T 



^n{LT)\ 
2,T j 



+ 



TTP,n 
-"2,T 



(17) 
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where rp is given by Eq. p5p with the LT term in the "quark-parton model" , that is with- 
out QCD radiative corrections to the coefficient functions. The deuteron structure functions 
calculated using Eq. ()17|) require 3-dimensional integration only. This approximation intro- 
duces only a marginal bias in the final results. 

We parametrize the functions H2^t{x) at the selected values of x and interpolate between 
those grid points using cubic splines. The positions of the grid points were selected in such 
a way to provide the overlap between the error-bars of nearest grid points. The values of the 
functions H2^t{x) at the grid points have been ffited to data. This method made it possible 
to describe the different structure functions in our ffi and, at the same time, to keep the 
number of fitting parameters reasonable. 

All experimental errors in the data including uncorrelated statistical, correlated system- 
atical, and errors in overall normalizations have been taken into account in the analysis 
using the covariance matrix approach. The error bands throughout the paper are due to 
linear propagation of the errors into the fitted parameters. In several of the experimental 
data sets used in our analysis the overall normalization has been fixed by the authors from 
comparisons to other data. 

For this reason we introduce an additional normalization parameter for each target (in 
some cases also for each energy) in the experiment and fit these parameters simultaneously 
with the parameters of PDFs and the HT terms. The overall normalization errors for such 
data subsets are accounted by the error propagation in the corresponding normalization 
factors through the general Hessian matrix of the ffi. Such treatment of the normalization 
errors allows for a maximal self-consistency of the analysis, incorporating all existing in- 
formation and minimizing the normalizaton uncertainties in the data. The re-normalized 
data subsets and their re-normalization factors, derived in the fit, are listed in Table HI The 
re-normalization factors are generally close to 1 for the deuteron data and somewhat higher 
for the SLAC proton data. 

B. Results 

The isospin asymmetries of the HT terms in F2^t obtained from this fit are shown in 
Figs. Eltoinias full lines, which correspond to the la bands of the total experimental errors 
|49|. In order to study the sensitivity of the results to various theoretical assumptions and 
estimate theoretical errors, we performed a number of fits using different approximations 
and assumptions. 

The central value results of these fits are also given in Figs. El to IHl Since the errors in 
i/2,T obtained in all variants of the ffi are correlated, the statistically significant shifts are 
those whose central values lie out of the bands given by the solid curves. 

For all variants of the fit we put the constraint H2,t{x) = as x ^ 0, since the analyzed 
data are not sensitive to the isospin asymmetry at small x. Relaxing this constraint does 
not improve the quality of the fit and in this case the values of if2,T(0) are comparable to 
zero within errors. 

The sensitivity of the isospin asymmetry of the HT terms to the treatment of Fermi 
Motion and nuclear Binding (FMB) in the deuteron is shown in Fig0 One observes that 
different approaches lead to variations in ~^ ^^P to several standard deviations at x > 0.75. 
The high sensitivity of the fit to details of the treatment of Fermi motion indicates that an 
accurate account of nuclear smearing in deuterium is crucial for the determination of i^^^^ 
at large x. 
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TMC also strongly affect the extraction of the isospin asymmetry due to their interplay 
with nuclear corrections at large x, as it can be seen from Fig|21 This is because TMC 
modify the x dependence of the LT structure functions at large x. The impact of TMC on 
the deuteron correction is illustrated in FigEJ in which the fit without TMC in pP'^^'^^^^ 
in Eq.(|17|) is compared to the fit with the full treatment of TMC. We observe a noticeable 



n—p 



difference between these two fits in the region of x ~ 0.75, especially for H2 

We have also studied the sensitivity of our result to the choice of the deuteron wave 
function in Eq.Q. Calculations of deuteron wave functions predict different amounts of high 
momentum components depending on the nucleon-nucleon potentials that are used, and on 
the treatment of relativistic effects. We considered two extreme situations by comparing in 
Fig|7| results using the Bonn ji^l and Paris i^l potentials that predict the smallest high 
momentum components among modern calculations, and a phenomenological wave function 
that reproduces y-scaling data and that has therefore a larger amount of high momentum 
components. Relativistic calculations seem also to have a larger amount of high momentum 
components ji^l. One can see that the functions ifg? vary within one standard deviation 
for the Paris and Bonn wave functions. The phenomenological distribution, however, is more 
than one standard deviation away at x > 0.75. 

The isospin asymmetry in the HT terms is also affected by off-shell effects in the bound 
nucleon structure functions discussed in Section ITTTl Since the calculation of off-shell effects 
is model dependent, these effects are in principle the main source of theoretical uncertainty 
in our analysis. However, the shift from the mass shell for the bound proton and neutron 
in the deuteron, measured by A = (p^ — M^)/M^, is small because of the weak binding 
in the deuteron (see Fig|2I). As a result, the net effect of off-shell corrections is within the 

? error-bars (see Fig|HJ and thus the uncertainty from the modeling of off-shell effect is 
effectively small for the deuteron. We note, however, that the off-shell effects will be much 
more important in heavy nuclei, as also shown in Fig|31 

In order to examine the impact of different models of nuclear effects in DIS, we repeated 
the fit with the deuteron correction calculated using the Nuclear Density Model (NDM) 
discussed in Section 11111 Since the NDM correction is not available for x > 0.8, we have 
removed the corresponding data points from the fit. Furthermore, since the NDM does not 
distinguish among structure functions, we assumed the correction to Ft to be the same as for 
F2. The results are shown in Fig|Sl The term HJ^~^ obtained from the NDM fit is different 
from the corresponding term of the FMB fit, while for H2~^ we observe a good agreement 
between these two fits. However, since the nuclear density model of the EMC effect in the 
deuteron is essentially different from the approach discussed in this paper, this agreement 
appears to be accidental. 

Another source of theoretical uncertainty comes from higher order QCD radiative correc- 
tions. These corrections decrease as increases. For this reason the radiative corrections 
can simulate the power-like terms in some kinematical regions of Q^. Indeed, it is well known 
that the magnitude of phenomenological HT terms is strongly correlated with the order of 
the QCD radiative corrections applied in the analysis. The HT terms drastically decrease 
when going from LO to NLO in the structure function fit. However, the magnitude of the 
HT terms does not change much when going from NLO to NNLO, the variation of the HT 
terms stays within one standard deviation [2, H^- In the present analysis we also observe 



only a marginal variation of ? after the NNLO corrections have been applied (see FigEl). 
We do not consider the soft gluon re-summation j2l| as well as the log-like dependence of 
H2T because of anomalous dimensions The isoscalar part of the HT terms can be 
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affected by these effects at large x ~ 0.9. However, the isovector combinations -ff^^^ seem 
to be much less sensitive to these effects. 

To summarize, the isospin asymmetry is stable with respect to the theoretical uncertain- 
ties addressed in this paper. The statistical correlations between the LT and the HT terms 
in our fit are illustrated in FiglTUl A wide kinematical region of considered DIS data allows 
to reliably separate the LT and HT terms and, as a result, the corresponding correlation 
coefficients are less than 0.5. 



V. PHENOMENOLOGY 

It is instructive to compare our results to the earlier extraction of the isotopic effects 
in the HT terms. The NMC extraction of the HT asymmetry in F2 [i^ is based on the 
equation 

Fl_ r,{F,"(^^)} / _ ci-cf] 

in which the combined SLAC-NMC-BCDMS data for the deuteron and proton are used. 
The relation between the functions Cf'" and the functions ifg'"" from our analysis can be 
written as follows 

ttP Tin 

It must be noted, that the correspondence between these two definitions is somewhat uncer- 
tain since the denominators in Eq. p8|) depend on Q^. However for comparison, in Figllllwe 
plot the difference — C2 calculated from our results using Eq. lfTHj) at different Q'^. Both 
extractions agree within errors. We note that our errors are systematically smaller because 
of a larger number of data points used in the analysis. It must be also emphasized that our 
analysis and the one by NMC are different in a few aspects (no deuteron corrections were 
applied in the NMC analysis, different treatments of systematic errors were used, etc.). All 
of these factors could be responsible for the remaining discrepancies. 
In Fig |12[ the isospin asymmetry in the ratio R = ol/ ctt 

R = 7'^^ - 1. 



extracted from the fit of the SLAC data |46[, is compared to our results. A good agreement 
on — BP suggests that our fit is self-consistent, since our analysis includes the data used in 
Ref. 113 . The value of - Rp at x = 0.03 ^ 0.35 measured by NMC ^ is also comparable 
to the calculations based on our fit. In Fig^Jwe show the different contributions to R^ — Rp 
considered in our analysis. It can be seen that at large x the value of R^ — R^ is defined 
mainly by the HT terms. 

The isotopic asymmetries i/g ? determined from our fit can be also compared with the 
predictions of different theoretical models. A popular model is the infrared renormalon 



model (IRR) of Refs. |lOl. In this model HT terms derive from the resummation of 

multi-loop diagrams, and their x dependence is obtained as the Mellin convolution of the 
LT terms with ffavor-independent coefficient functions C^^^, 



13 



The dimensional normalization factor A2 determines the characteristic scale of the HT terms. 
This scale is not determined in the IRR model and it is adjusted from comparisons with the 
data. In Fig lT^ our results are compared to the IRR model calculation for the non-singlet 
contribution, with the normalization factor A2 = —0.3 GeV^. It is clear that the agreement 
is poor and it cannot be improved by simply varying the normalization factor of the IRR 
model. Alternatively, the comparison can be made in terms of the Mellin moments 

Jo 

The Mellin moments of -ff^^^ are given in Table IVII In Fig lT^ the ratio of these moments 
to the moments of p^^^ is plotted. Also shown are the moments of C^^^. We only 
show N > 2, since the data do not allow us to constrain the behavior of small x. 

Again, we observe a disagreement with the IRR model. Par ameterizat ions of DIS structure 
functions interpolating between the low and high regions are also available 0, 0] . The 
comparison of our results with these models is not directly feasible, because of the strong 
dependence introduced in their evaluation of the HT terms, due to effects beyond the 
OPE. 

Having estimated the power corrections to different structure functions, one can make 
an extrapolation into the resonance region in order to study the phenomenon of duality 
(48| . Examples of such extrapolations are given in Figs llSllTl At values of relevant for 
the experiments at Jefferson Lab the HT terms contribute moderately to the ratio F2/F2. 
For the difference — the impact of the HT terms is larger. At < 1.4 GeV the 
dominant contribution comes from ifg"^. At > 1.4 GeV the contribution from ifg"^ is 
small and the main effect comes from HJ^~^. However the errors are large in this region, 
since it corresponds to a; ~ 0.4, where HJ^~^ has its largest uncertainty. 

VI. CONCLUSIONS 

In conclusion, we determined the isospin asymmetries of the functions H2 and de- 
scribing the HT terms of the DIS structure functions. The value of is consistent with 
zero within the errors for all values of x. Also H2 is consistent with zero at low and inter- 
mediate values of x; It deviates from zero at a; > 0.7. We performed a careful study of the 
theoretical uncertainties that might affect the extraction and we conclude that they do not 
overwhelm the effect. The asymmetry H2~^ is negative at large x. It reaches its maximum 
at X ~ 0.8 where it is ~ 0.03 GeV^, in agreement with the order of magnitude of the scale 
of QCD, A^. We also find that the x dependence of -ff2T^ poor agreement with the 
predictions of the IRR model. For more conclusive comparisons more precise data at a; ~ 0.4 
and ~ 1 GeV^ are necessary. 
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TABLE I: List of 


re-normalized experiments with the corresponding 


; re-normalization factors r/. 


Experiment 




r?[%] 




proton 


deuterium 


SLAC-E-49A 


1.5 ± 1.3 


-1.1 ± 1.2 


SLAC-E-49B 


2.8 ± 1.3 


0.2 ± 1.2 


SLAC-E-87 


2.9 ± 1.2 


0.9 ±1.2 


SLAC-E-89B 


1.4 ± 1.2 


-1.1 ± 1.2 


SLAC-E-139 




0.6 ± 1.3 


NMC(90 GeV) 


-0.8 ± 1.4 


-2.3 zb 1.3 


NMC(120 GeV) 


0.8 ± 1.3 


-1.4 ± 1.3 


NMC(200 GeV) 


2.4 ± 1.3 


0.1 ± 1.3 


NMC(280 GeV) 


1.2 ± 1.3 


-1.0 ± 1.3 


TABLE IL The Mellin moments of 


N 




M^-P(iV) 


2 


-0.0058 ± 0.0069 


-0.012 ±0.014 


3 


-0.0046 ± 0.0024 


-0.0048 ± 0.0052 


4 


-0.0041 ± 0.0013 


-0.0020 ± 0.0026 


5 


-0.00355 ± 0.00084 


-0.0008 ±0.0015 


6 


-0.00301 ± 0.00061 


-0.00019 ± 0.00096 


7 


-0.00254 ± 0.00048 


0.00009 ± 0.00068 


8 


-0.00215 ± 0.00040 


0.00023 ± 0.00052 


9 


-0.00183 ± 0.00034 


0.00029 ± 0.00041 


10 


-0.00157 ± 0.00030 


0.00030 ± 0.00034 
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FIG. 1: Deuteron Compton scattering amplitude in incoherent scattering approximation. 
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FIG. 2: The ratio R2 calculated in different approximations. In the upper panel this ratio is 
presented as a function of x for fixed Q^: Fermi motion and binding effects (dotted line), Fermi 
motion and binding effects combined with target mass corrections (dashed line) ; the full calculation 
including Fermi motion, binding, target mass and off-shell corrections is given by the solid line. 
The shaded area in the upper panel corresponds to the prediction of the nuclear density model of 
Ref . [3^ . In the lower panel the ratio R2 is shown as a function of W for a few different . 
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FIG. 3: Average nucleon kinetic and separation energy and off-shellness A as functions of the 
nuclear mass number A. 
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FIG. 4: Kinematic region of the data used in analysis. The curves correspond to constant values 
of the invariant mass W whose values in units of GeV are indicated in the plot. 
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FIG. 5: Isospin asymmetries of the HT terms obtained using different treatments of Fermi motion 
and binding corrections: Eqs.Q (delimited by solid lines, see text); Atwood-West ji^ (dash- 
dotted line); Eqs.@ in the Bjorken limit, i.e. if all l/Q^ terms were disregarded (short dashes). 
The curve with long dashes shows the result without Fermi motion and binding corrections. 
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0.04 




FIG. 6: Sensitivity of the isospin asymmetries to different approximations used in the calculation 
of the deuteron correction. Dashed line: fit with no TMC in the calculation of the deuteron 
correction; Area between solid lines: fit with full treatment of the TMC). 
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FIG. 7: Results of the fits with different deuteron wave functions: Ref.j^^] (area between solid 
lines); Ref. 41 1 (dashed line); Ref.ji^ (dashed-dotted line). 
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FIG. 8: Isospin asymmetries in the HT terms obtained within the incoherent scattering approxima- 
tion (area between sohd hnes) as compared to the results obtained with the treatment of ofF-shell 
corrections (dashed-dotted hnes). Resuhs of the analysis based on the nuclear density model of 
Ref.js^ are given by the dashed lines. 
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FIG. 9: Impact of NNLO QCD corrections on the isospin asymmetries of the HT terms: NLO 
QCD fit (area between sohd hnes); NNLO QCD fit (dashed hnes). 
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FIG. 11: The isospin asymmetry of the HT coefficients determined by NMC (points with errors 
bars) compared to to our determination (area inside error bands). 
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FIG. 12: The isospin asymmetry in the structure function R determined from the SLAG data 
(points with error bars) compared to to our results. The LT contribution corrected for the TMG 
with no deuteron correction (dotted hue); both the LT and HT terms but no deuteron corrections 
(dashed hne); the same as above with the deuteron correction (area between solid lines). 
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FIG. 13: Comparison of the phenomenological isospin asymmetries (area between solid lines) to 
the prediction of the infrared renormalon model (dashed lines). 
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FIG. 14: The ratios of the MeHin moments of HT terms to the ones of LT terms (points with 
error bars) are compared with predictions of the IRR model (area between soHd hnes) . The inner 
bars give total experimental errors in the moments, the outer bars include the error due to the 
extrapolation into the unmeasured region at x > 0.9. 
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FIG. 15: The ratio jF^ calculated in different approximations and extrapolated to the resonance 
region: LT terms only (dotted line); effect of TMC (dashed line); TMC and HT terms (area between 
solid lines). 
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FIG. 16: The difference R — BP calculated in different approximations and extrapolated to the 
resonance region: only the LT terms with the TMC (dotted line); the same as above and the 
contribution from H2 (dashed line); the same as the dashed line and the contribution from i^x 
(area between solid lines). 



33 



